Gravitational waves from color-magnetic 'mountains' in neutron stars 
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Neutron stars may harbour the true ground state of matter in the form of strange quark matter. 
If present, this type of matter is expected to be a color superconductor, a consequence of quark 
pairing with respect to the color/flavor degrees of freedom. The stellar magnetic field threading 
the quark core becomes a color-magnetic admixture and, in the event that superconductivity is 
of type II, leads to the formation of color-magnetic vortices. In this Letter we show that the 
volume-averaged color-magnetic vortex tension force should naturally lead to a significant degree 
of non-axisymmetry in systems like radio pulsars. We show that gravitational radiation from such 
color-magnetic 'mountains' in young pulsars like the Crab and Vela could be observable by the 
future Einstein Telescope, thus becoming a probe of paired quark matter in neutron stars. The 
detectability threshold can be pushed up toward the sensitivity level of Advanced LIGO if we invoke 
an interior magnetic field about a factor ten stronger than the surface polar field. 
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Introduction. — Neutron stars might actually be quark 
stars in the event that stable strange quark matter is 
realised in the high density / pressure environment of these 
objects [H. Confirming (or ruling out) the existence of 
this type of matter ranks among the ultimate goals of 
neutron star astrophysics. But this will not easily come 
about; quark stars may masquerade as neutron stars with 
respect to the bulk properties of mass and radius [3] . The 
situation is exacerbated if the star is a hybrid, with a 
uark inner core cloaked by a mantle of hadronic matter 
Probing, for instance, the physics of the crust may 
reveal nothing about a quark core. Similarly, if stellar 
pulsation modes were to be observed, their properties 
could be controlled by the hadronic outer core and crust 
e.g. this is likely the case for the r-mode instability [4, 
. Neutron star cooling, a potentially powerful probe of 
stellar structure, could also prove to be largely insensitive 
to the physics of an inner quark core @ . 

If present, quark matter is likely to be in its stablest 
form as a condensate of paired quarks. Among the var- 
ious possible states, the two most intensively studied 
are [6|: (i) the color-flavor- locked (CFL) phase where all 
three quark species {u, d, s) are paired (ii) the two-flavor 
phase (2SC) where only the u and d quarks pair. In both 
phases the system is a 'color superconductor', an exotic 
variant of ordinary electric superconductivity. 

This key property of quark matter is the subject of 
the present Letter. The large scale magnetic field present 
in neutron stars can penetrate the CFL/2SC region, in- 
ducing 'color-magnetic' stresses which deform the stel- 
lar core, much like the magnetic deformation in ordinary 
neutron stars [31 ■ As we discuss in detail below, in the 
event that the system behaves as a type II superconduc- 
tor the deformation is greatly amplified by the formation 
of color-magnetic vortices, each one carrying an energy 



^ 10"^ times higher than the energy of ordinary protonic 
vortices (a result of the stronger Meissner screening as- 
sociated with color superconductivity). 

If systems like radio pulsars do contain a color- 
superconducting quark core, threaded by a large-scale 
magnetic field misaligned with the rotation axis, they 
should display a significant non-axisymmetric deforma- 
tion in the form of a color- magnetic 'mountain'. As a 
result of this, these systems will emit a continuous grav- 
itational wave signal, at a level much stronger than that 
due to magnetic field deformation in their hadronic coun- 
terparts. Our results on the detectability of such signals 
suggest that young pulsars - most notably the Crab - 
could be observed by future gravitational wave detectors 
like the Einstein Telescope (ET) and possibly even Ad- 
vanced LIGO. 

Color superconductivity in neutron stars. — A mag- 
netic field threading a region of color-superconducting 
quark matter is 'rotated' in the space of internal gauge 
degrees of freedom, and becomes an admixture of 'color- 
magnetic' fields (see Ref. [8] for a review). Among the dif- 
ferent possible combinations, the following mixed gauge 
field A"^ is associated with a massive propagator and su- 
perconducting Meissner screening [8] 

= sinxA + cosxA^ (1) 

where A and A* are the electromagnetic and (one of) the 
gluonic vector potentials, respectively. The 'mixing' an- 
gle X is (essentially) the electromagnetic to gluonic cou- 
pling constant ratio (expressed in ?i = c = 1 units) 
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(l + 3gV4e2)-i/2 (CFL phase), 
(l + 3<7Ve2)-i/2 (2SC phase). 



(2) 



The slightly different numerical factors in these expres- 
sions refiect the different quark pairing in the two phases. 
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The disproportionate coupling strengths, « 47r/137 
and g « 3.5 in natural units, result in % <C 1 which 
in turn implies that the field is mostly gluonic. The 
bulk electromagnetic field is insensitive to the quark color 
pairing/superconductivity and the system behaves para- 
magnetically with respect to it [6, 8]. 

The potential ([1]) is associated with a color-magnetic 
field B-^ = V X A-^ which obeys an Ampere-type law 
V X B''^ — Ati3^/c where J-^ is the combined color- 
electric current i9| . In the context of laboratory superflu- 
idity/superconductivity, the mass/electric current makes 
direct contact with the macroscopic quantum state of 
the paired system. The same is true for the more ex- 
otic quark color superconductor. The current has a form 
J-^ = ki\79 + k2A^ where fci, ^2, depend on the super- 
conducting phase (see for the explicit expressions). 
The (spatially varying) 9 is the so-called order parameter, 
that is, the phase of the macroscopic wavefunction. 

Color superconductivity displays the familiar type I/II 
dichotomy. Which type prevails in the conditions met 
in neutron stars is largely determined by the pairing gap 
energy to quark chemical potential ratio A/ fiq [TTl . The 
formal criterion for determining the type of superconduc- 
tivity is the same as in electromagnetic systems: type II is 
the energetically favorable state provided the Ginzburg- 
Landau parameter k = A/^ > l/\/2, where A and ^ are, 
respectively, the penetration and coherence lengths. The 
situation is blurred by the fact that the underlying theory 
is actually an approximation near the critical tempera- 
ture Tc for the onset of superconductivity. Extrapolating 
to the temperature regime T Tc (appropriate for as- 
trophysical neutron stars) for the 2SC phase, Alford & 
Sedrakian 11| estimate k w llA//iq which predicts a 



type II state for A > 25 MeV (assuming a typical quark 
chemical potential fiq = 400 MeV). The theoretically 
predicted gaps lie in the range A ^ 20 — 100 MeV [6| and 
therefore most of them satisfy this condition. The anal- 
ysis for the CFL state is equally uncertain. The results 
of [I4I suggest K sa 9A//^q which again predicts type II 
superconductivity for 'typical' pairing gaps. 

Without any additional microphysics input we here- 
after assume that CFL/2SC color superconductivity is of 
type II and that the B-^ field forms a stable Abrikosov 
lattice of color-magnetic vortices. The formation of vor- 
tices manifests in the form of singularities in the order 
parameter 9 at the location of each vortex. Skipping the 
details that are not essential for our present discussion, 
the end-product is the formulation of the following quan- 
tization condition for an isolated vortex (see jlo|): 
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where the asterisk stands for ★ — {2SC, CFL} and 0^ is 
the color-magnetic flux quantum (see below). The vor- 
tex singularity is located at r = with the vortex itself 
identified with z-axis in cylindrical coordinates {r, ip, z}. 



The condition (j3]) applies on the mesoscopic neighbor- 
hood of a single vortex, where the matter can be assumed 
to be uniform (this entails A^, « constant). Taking the 
curl of ^ leads to the London equation 



V^B^ - -Vb^ 

A? 



(4) 



This equation encapsulates the Meissner-screening prop- 
erty of the B^ field. 

The color-magnetic vortex energy. — Using the im- 
posed symmetries of our single vortex system the above 
London equation reduces to an ODE with the familiar 
Bessel-function solution Il3t. 
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This field is associated with the flux J dS -B^ — (/)^. The 
total energy per unit vortex length consists of contribu- 
tions from the color-magnetic field and current. 
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It is straightforward to compute the integrals using 
If we furthermore assume a strong type II superconduc- 
tivity, i.e. <C Aj;, the final result is 
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Not surprisingly, this energy has a functional form iden- 
tical to that of fluxtubes in a proton superconductor [14 1 . 

We now need to input numerical values for A^, and (f)^,. 
These are given in [10| ; the penetration length is 

A.«lG-i3^,(l/50/i^lioo cm (8) 

where /iq4oo = ^q/400MeV, Ccfl = 1-1, C2SC = 1-9, 
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52SC = o + e 



, 5CFL = - (3g + 4e ) 



(9) 

Note that these coupling combinations are to be evalu- 
ated using natural units. For the flux quantums we have 



4>i, = hc/g^, = 6smx(t>o, 



(10) 



where (po = hc/2e w 2 x 10-"^ G/cm^. 

Using these estimates and making the usual approxi- 
mation (1/2) log(A^,/^i,) « 1 [13], we arrive at our final 
results for the color-magnetic vortex energy: 



4.9 X 10-*/Zq4oo ergs/cm. 



£CFL _ 1 c inl0,,2 



1.5 X 10^'^/iq4no ergs/cm. 



(11) 



(12) 



Most noteworthy, the vortices of both CFL and 2SC 
phases carry much more energy than the energy £p w 
10^ ergs/cm of protonic fiuxtubes [l4| . 
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The stellar deformation. — We can estimate the elliptic- 
ity due to the color-magnetic stresses in the quark core 
using the ratio of vortex to gravitational binding energy: 



ex 



3GAP/{4R)' 



(13) 



where J\fx is the vortex area density averaged over the 
quark core, Vq is the volume of the core, and we have 
used the gravitational binding energy appropriate to an 
n = 1 polytrope [hereafter, we scale our results to a mass 
M = IAMq, a radius R = 11.7 km, and a moment of in- 
ertia I = 10*^ g cm^, consistent for an n = 1 polytrope]. 
It should be pointed out that ([T^ is an approximation [7| ; 
nevertheless, it is known to be quite accurate in repro- 
ducing rigorously computed ellip ticities of magnetically 
deformed neutron stars (see [15| for a discussion of the 
superconducting case). 

Following Ref. [ll[, we can relate Afx to the macro- 
scopic interior magnetic field B: 

A/k = sinxB/(/i* = Afp/6 « 8.1 x lO^Bia cm'^ (14) 

where Afp — B / (pQ \s the density of protonic vortices and 
Bi2 = B/10^^ G. Parametrizing the quark core as = 
fVsta.1-, we obtain 
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FIG. 1: Estimated gravitational wave amplitudes for known 
pulsars, assuming they contain 2SC or CFL cores with color- 
magnetic ellipticities given by (fT5)) and (fT5)l . We assumed a 
quark volume fraction / — 1/2, a quark chemical potential 
of /iq = 400 MeV, and a ratio of internal volume-averaged to 
polar field of a = B / Bp = 2. For the Crab and Vela, the short 
vertical line segments show the range in wave amplitude as fj,q 
varies over the interval 300-600 MeV, assuming CFL cores. 
All pulsar data was taken from the ATNF database [IJ] . The 
detector noisecurves were taken from [2^ and [23| . 
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where B12 is the volume-average of B over the quark 
core. The CFL-phase ellipticity is similarly found to be: 



e^FL « 2.4 X 10-' /Bi2M^400^11.7^A" 



(16) 



For comparison, the ellipticity resulting from an ordinary 
protonic superconducting star is (setting / = 1), 



9.7 X IQ-^^BuRti.rM^l 



(17) 



Gravitational wave detectability. — It is straightforward 
to use these results to estimate the level of gravitational 
wave emission from known pulsars, assuming they do in- 
deed contain color superconducting 2SC or CFL cores. 

First we need to use the observed spin and spin-down 
rate (l to infer the magnetic field strength. Note that 
the volume-averaged field B is likely to exceed the dipole 
surface polar field Bp by a factor of a few; certainly this is 
the case for known magnetic equilibria, see e.g. 
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where B w 2Bp. We therefore define B = aBp. For the 



electromagnetic spin-down torque we will use the usual 
vacuum dipole resuU: Jew = {R'^ / 6c^) B^n^ = AB^n^. 
For the gravitational wave spin-down, we have Jqw — 



{32GP/5c^ 
and (H 



lein^ nj, which. 



takes the form Jqw — Cxa B^fl'>. 



when combined with (|15|) 
2n2o5 Then the 



balance Jew + >^gw = leads to 



-,1/2 



(18) 



Note that Jew in ([T5|) is negligible in all cases apart from 
rapidly spinning systems with a ^ 100. 

The field B = aBp can then be substituted into equa- 
tions (fT5)) . (fT6|) to give the corresponding ellipticities, 
which can then be converted into gravitational wave am- 
plitudes using equation (55) of Ref. [20| : 
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1/2 



(19) 



for a star a distance d from Earth. 

The results are shown in Figure [U for both 2SC and 
CFL cores. We have used sensible values for the vari- 
ous physical parameters, including a quark chemical po- 
tential — 400 MeV, a fiducial quark core volume 
fraction / = 1/2, an internal average magnetic field of 
a = B / Bp ^ 2. Note that the calculated ellipticities, and 
therefore the gravitational wave amplitudes, are propor- 
tional to affj,^, so that the results may scale up or down, 
depending upon the actual values of these parameters as 
realised in Nature. To capture some of this uncertainty, 
for the CFL estimates for the Crab and Vela pulsars, we 
have held all parameters fixed apart from /^q, which we 
have varied over the interval 300-600 MeV [6|, to give 
the two short line segments indicated in the Figure. 

Discussion. — Clearly, the Crab and Vela pulsars are 
the two most interesting candidates in terms of detection; 
both lie close to the ET detectability curve. In the case 
of the Crab and a CFL core, our results suggest that, 
for the stellar parameters used here, the Crab may be 
detectable by ET. In the case of a 2SC core, detectability 
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TABLE I: Minimum value of a = B /Bp, the ratio of volume- 
averaged internal field to polar field, required for gravitational 
wave detection. We fix / = 1/2 and = 400 MeV. 





ET 


Advanced LIGO 




Crab 


Vela 


Crab 


Vela 


CFL 


0.78 


3.3 


12 


46 


2SC 


2.3 


9.9 


35 


137 



may be possible, but only if /^q lies in the upper part of 
the range quoted above. In the case of the Vela and 
a CFL core, detection would be possible only if /iq lies 
close to the upper end of the quoted /iq-range, while a 
2SC core in Vela would not be detectable. Clearly, the 
systems lie sufficiently close to the detectability curve 
that a confident statement regarding detectability cannot 
be made here; a more careful treatment is clearly needed. 

There has been discussion as to the possibility of some 
pulsars harbouring much stronger fields than the spin- 
down estimates would suggest. In particular, the ex- 
istence of 'anomalous' braking indices ClQ./CP' < 3 has 
been taken as evidence of the gradual emergence of a 
buried strong magnetic field from the interior of some 
young pulsars by diffusive processes (see 12J| and refer- 
ences therein). In fact, both Crab and Vela belong to 
this pulsar subgroup. Further weight has been added to 
this theory by observations of the rotation-powered pul- 
sar PSR J1734, which seems to be evolving toward the 
magnetar portion of the pulsar (r2,f2) plane [IJ. This 
would correspond to the young pulsars having internal 
fields one to two orders of magnitude larger than their 
polar fields, i.e. a ^ 10-100. 

In light of these suggestions, we give in Table |T] the 
minimum a values required for detectability of the Crab 
and Vela in our model, assuming sensible values of all 
other parameters (/ = 1/2, ^q = 400 MeV). Clearly, if 
the Crab and Vela do harbour strong internal fields, the 
corresponding CFL and 2SC emission would lie within 
range of ET, and possibly even Advanced LIGO. We also 
note that the current 95% confidence direct upper limit 
on gravitational wave emission from the Crab pulsar of 
/j95% ^ 2.0 X 10-25, obtained using LIGO data Tt^, im- 
plies an upper limit of a ~ 250 for this pulsar, assuming 
that it contains a CFL core. 

In fact, for a ~ 100, we find that ~ 40 pulsars are 
detectable by ET, including many millisecond pulsars. 
Increasing a to ^ lO'^ pushes many of the millisecond sys- 
tems close to their spin-down limits, leaving ~ 30 of them 
detectable by ET (such high a values could naturally 
materialise in recycled millisecond pulsars as a result of 
magnetic burial during the previous accretion-dominated 
stage of their lives). For these systems, the gravitational 
wave torque would then dominate the electromagnetic 
one, and the conventional electromagnetic-only calcula- 
tion of Bp would significantly overestimate the actual 



polar field strength (see equation (IT8l) l. 

Failure of ET to detect gravitational waves from the 
Crab or Vela does not rule out the existence of color 
superconducting quark cores in these (or any other) neu- 
tron stars. It may be that the relevant parameters lie at 
the low end of the possible ranges, i.e. /J,q ss 300 MeV, 
and/or the core occupies only a small fraction of the total 
stellar volume. Alternatively, it may be that the internal 
field arranges itself in such a way as to give a smaller ellip- 
ticity than our estimate of equation ([T3|) would suggest. 
However, this would require a rather special geometry, 
not seen in magnetic fields studied in 'normal' neutron 
stars (see e.g. [l8j). Another possibility is that the vortex 
arrays discussed here simply aren't stable, in other words 
type II color-superconductivity is not realised, resulting 
in a lower ellipticity than equation (|13|) would predict. 



However, if the link between field strength and defor- 
mation is similar to the well understood problem in non- 
color superconducting stars, our estimates should be re- 
liable. In this sense the link between CFL/2SC cores 
and gravitational wave emission is relatively robust, and 
more direct than the link between exotic elastic cores 
and gravitational wave emission, where, in addition to 
the existence of the exotic core, a mechanism to generate 
the non-axisymmetric strains that perturb the star away 
from the lower energy axisymmetric state must also be 
provided 



A more accurate calculation of the gravitational 
wave emission would self-consistently solve for a self- 
gravitating color superconducting star, accounting for 
the multi-component nature of the stellar matter and the 
magnetic field geometry, and differentiating between the 
quark core and hadronic envelope (see [26j for discus- 
sion of this last point). Equally crucial is the need for a 
more accurate input from microphysics, e.g. vortex struc- 
ture calculations for the density /temperature regime ap- 
propriate to neutron stars and consideration of a larger 
variety of color superconducting states While chal- 
lenging, the proximity of candidate sources to the ET 
detectability curve clearly motivates the development of 
the tools needed to do this, to better assess the ability 
of gravitational wave observations to provide this unique 
probe of the state of matter at the highest densities. 
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